function fval = average_log_normal(s, c_mu, c_sig)

%int_method = 'auto';
Rel_Tol    = 1e-5;
Abs_Tol    = 1e-6;
fval = integral(@(x)1 / (c_sig * sqrt(2 * pi)) * exp(- (log(x) - c_mu).^2 / 2 / c_sig^2), 0, s, 'RelTol', Rel_Tol, 'AbsTol', Abs_Tol);

end